How Perfect a Gluon Plasma Can Be in Perturbative QCD? 
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The shear viscosity to entropy density ratio, rj/s, characterizes how perfect a fluid is. We calculate 
the leading order rj/s of a gluon plasma in perturbation using the kinetic theory. The leading order 
contribution only involves the elastic gg gg (22) process and the inelastic gg ggg (23) process. 
The Hard- Thermal-Loop (HTL) treatment is used for the 22 matrix element, while the exact matrix 
element in vacuum is supplemented by the gluon Debye mass insertion for the 23 process. Also, 
l" ■ the asymptotic mass is used for the external gluons in the kinetic theory. The errors from not 

implementing HTL and the Landau-Pomeranchuk-Migdal effect in the 23 process, and from the 
^ . uncalculated higher order corrections, are estimated. Our result for rj/s lies between that of Arnold, 

■ Moore and Yaffe (AMY) and Xu and Greiner (XG). Our result shows that although the finite 
angle contributions are important at intermediate a s (ct s ~ 0.01-0.1), the 22 process is still more 
important than 23 when a s < 0.1. This is in qualitative agreement with AMY's result. We find no 

K-j . indication that the proposed perfect fluid limit rj/s ~ l/(47r) can be achieved by perturbative QCD 

■ alone. 

43 I. INTRODUCTION 

PK 

A perfect fluid is a system with zero shear and bulk viscosities, rj and £, and no dissipation. These conditions can be 
satisfied for a superfluid at zero temperature where only the super fluid component exists, but a sharper description 
is with the dimensionless ratios rj/s and (/s, where s, the entropy density, vanishes for the superfluid component as 
| well. While £/s can still be zero for scaling invariant systems, the situation for rj/s is more subtle, 
i In general, stronger interactions implies a smaller r\. Thus, a perfect fluid with the smallest rj/s is likely to be 
strongly interacting which requires nonHDerturbative tools to compute it. The anti-de-Sitter space/conformal field 
theory correspondence (AdS/CFT) [l|-[3| allows the rj/s of strongly interacting CFT's to be computed in weakly 
interaction gravitational theories. A universal number rj/s = l/(47r) is found for every CFT with a gravity dual in 
the large N, with N the size of the gauge group, and infinite t'Hooft coupling limit [J-Q- With this result, together 
• ■ with the connection to the uncertainty principle through the relation r// s ~ AEAt, with AE and At the mean energy 
and life time of quasiparticles, Kovtun, Son, and Starinets (KSS) Q conjectured that the strongly interacting CFT 
value l/(4.7r) is the minimum bound for rj/s for all physical systems. 

Theoretically, there are several attempts to evade this bound. It is found that rj/s can be as small as possible (but 
• • . still non- negative) in a carefully engineered meson system [jj, |8|, although the system is meta-stable. Also, in strongly 
interacting CFTs, l/N corrections can be negative [H [l(| and can modify the rj/s bound slightly [Tit fT3| . 

Experimentally, there are intensive interests to find the most perfect fluid (see [H, [H| for recent reviews). The 
<_i smallest rj/s known so far is realized in a system of hot and dense matter thought to be quark gluon plasma just 
above the phase transition temperature produced at RHIC ^5T4l7| with rj/s = 0.1 ± O.l(theory) ± 0.08(experiment) 
[l8| . A robust upper limit rj/s < 5 x l/(47r) was extracted by another group [l!| and a lattice computation of gluon 
plasma yields rj/s = 0.134(33) [20| . Progress has been made in cold unitary fermi gases as well. An analysis of the 
damping of collective oscillations gives rj/s > 0.5 [2ll[2^. Even smaller values of rj/s are indicated by recent data on 
the expansion of rotating clouds [23|, [24j but more careful analyses are needed [25ll2r|. 

Even if the l/(47r) bound for rj/s turns out to be invalid, it is still interesting to use it as a bench mark value for 
the perfectness of fluids. It was found that based on the perturbative QCD (PQCD) analysis of Arnold, Moore and 
Yaffe (AMY) [27],[28j], the measured rj/s at RHIC cannot be explained by PQCD (for a recent review, see, e.g., [29|). 
This strongly interacting QGP picture is very different from the conventional picture of weakly interacting QGP and 
is considered as one of the most surprising discoveries at the RHIC. 

However, a recent perturbative QCD calculation of rj/s of a gluon plasma by Xu and Greiner (XG) [3(| shows 
that the dominant contribution comes from the inelastic gg ■<-> ggg (23) process instead of the elastic gg —¥ gg (22) 
process. In particular, the 23 process is 7 times more important than 22. Thus, rj/s ~ 1/47T can be achieved when 
the strong coupling constant a s ~ 0.6. In Ref. [3l|, XG and their collaborators improve their calculation using the 
Kubo relation and give smaller contribution from the 23 process: about 5 times (2^9 times) the 22 process. Thus, 
the conventional weakly interacting QGP could still be valid. This is in sharp contrast to AMY's result where the 23 
process only gives ~ 10% correction to the 22 process. 
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Both XG and AMY use kinetic theory for their calculations. The main differences are (i) XG uses a parton cascade 
model |32| to solve the Boltzmann equation and, for technical reasons, gluons are treated as a classical gas instead of 
a bosonic gas. On the other hand, AMY solves the Boltzmann equation for a bosonic gas. (ii) AMY approximates the 
Ng f-> (N + l)g processes, N = 2, 3, 4 . . ., by the g «-» gg splitting in the collinear limit where the two gluon splitting 
angle is higher order. XG uses the soft gluon bremsstrahlung limit where one of the gluon momenta in the final state 
of gg —> ggg is soft but it can have a large splitting angle with its mother gluon. 

In an earlier attempt to resolve the discrepancy between XG's and AMY's results [33|, a Boltzmann equation 
computation of r/ is carried out without taking the classical gluon approximation (like AMY's approach) but the 
soft gluon bremsstrahlung limit is applied to the 23 matrix element (like XG's approach, modulo a factor 2 in the 
23 matrix element squared; sec 33] for details). It was found that the classical gas approximation does not cause 
a significant error in r\js (although the individual errors on rj and s are larger). However, the result is sensitive to 
whether the soft gluon bremsstrahlung limit is imposed on the phase space or not. If this limit is imposed, the result 
is closer to AMY's; if not, the result is closer to XG's. This raises the concern whether this approximation is good 
for computing rj. 

The goal of this paper is to settle this issue by removing both the soft gluon bremsstrahlung approximation and 
the collinear approximation to the 23 process. The leading order [0(aj 2 )] contribution to 7? only involves the 22 
and 23 processes [28[ (the power counting for 22, 23 and other processes are reproduced in 33j). In this paper, the 
Hard- Thermal-Loop (HTL) treatment is used for the 22 matrix element, while the exact matrix element in vacuum 
is supplemented by the gluon Debye mass insertion for the 23 process. Also, the Debye mass is used for the external 
gluon mass in the kinetic theory as well. The errors from not implementing HTL and the Landau-Pomeranchuk-Migdal 
effect in the 23 process, and from the uncalculated higher order corrections, are also estimated. 



II. KINETIC THEORY BEYOND THE SOFT OR COLLINEAR GLUON APPROXIMATIONS 

Using the Kubo formula, 77 can be calculated through the linearized response function of a thermal equilibrium 
state 
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At' I dt I d£ 3 ([T^(0),T^'(x,i)]), 
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where T lJ is the spatial part of the off-diagonal energy momentum tensor. In a leading order (LO) expansion of the 
coupling constant, there are an infinite number of diagrams [3~H . l35l. However, it is proven that the summation of the 
LO diagrams in a weakly coupled </> 4 theory [34 38] or in hot QED [39| is equivalent to solving the linearized Boltzmann 
equation with temperature-dependent particle masses and scattering amplitudes. The conclusion is expected to hold 
in weakly coupled systems and can as well be used to compute the LO trans por t coefficients in QCD-like theories 
[13, HI], hadronic gases [40l - |45| and weakly coupled scalar field theories (MSlH, S3 • 

The Boltzmann equation of a hot gluon plasma describes the evolution of the color and spin averaged gluon 
distribution function f p (x) which is a function of space-time x = (i, x) and momentum p = (E p , p). The infinitesimal 
deviation of f p (x) from its equilibrium value f^ q — (e v ' p ^ T — 1) _1 is denoted as 



X P (1 



where \p = xi x iP) can be parametrized as 
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at the leading order of the derivative expansion of the fluid velocity v(x) = (v , v). T = T(x) is the local temperature, 
p is the unit vector in the p direction. Bij(p) = B(p)(p l p : > — MjP) and B{p) are functions of p which will be 

fixed by the Boltzmann equation corresponding to the bulk and shear viscosities, respectively. In this work we will 
just focus on the shear viscosity calculation. 

The Boltzmann equation j48l453| for the gluon plasma reads 
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where the collision rates are given by 
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iV s = 16 is the color and spin degeneracy of a gluon. The i-th gluon is labeled as i while the ra-th gluon is labeled 
as p. For a process with I initial and (n — I) final gluons, the symmetry factor N(n, I) = l\(n — I — 1)!. For example, 
processes 12 ->■ 3p, 12 -s- 34p, 123 -J- 4p yield (n, Z) = (4, 2), (5, 2), (5, 3) and iV(n, = 2, 4, 6, respectively. 
In vacuum, the matrix element squared for the 22 process is 



|M 12 ^ 3 4| 2 = ^(4 7 r) 2 iV^l3-^- 
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where ot s = g 2 /(47r) is the strong coupling constant, and (s,t,u) are the Mandelstam variables s = (pi +P2) 2 , 
t = (pi - p 3 ) 2 and u = (pi - pi) 2 . 

For the 23 process [H], HU, under the convention Yli=iPi = 0, we have 



|-^12345-s-o| 2 — 1-^^0^12345 1 2 
^3 Ar2„ ,3 



= 54tt j 7V>; 



(12) 4 + (13) 4 + (14) 4 + (15) 4 



+ (24) 4 + (25) 4 + (34) 4 + (35) 4 + (45) 4 
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(23) 4 
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where (ij) = Pi-pj and the sum is over all permutations of {1, 2, 3, 4}. To convert to the convention P1+P2 — P3+P4+P5, 
we just perform the replacement: 
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In the medium, the gluon thermal mass effect serves as the infrared (IR) cut-off to regularize IR sensitive observables. 
The most singular part of Eq. (|6j) comes from the collinear region (i.e. either t ss or u ~ 0) which can be regularized 
by the Hard-Thermal-Loop (HTL) corrections to the gluon propagators [IB, E3] and yields [HH , 



M 12 ^ 34 | 2 « 7(12™ S ) 2 7V 2 (4£; 1 ^ 2 ) 2 
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where q = pi — p 3 = (qo,q), x = qo/\q\ and (j) is tne angle between pi x q and p 2 x q. The HTL self-energies 
(longitudinal) and (transverse) are given by 



Ht — mjj 



X 1 + X 7T_ 



y + fCl-^ln 1 



(10) 



The external gluon mass (i-e. the asymptotic mass) is the mass for an on-shell transverse gluon, and m 2 ^ — 
TIt {\x\ = 1) = mjj/2 both in the HTL approximation and in the full one-loop result. 

Previous investigations of the thermodynamics within resummed perturbation theory showed that the most impor- 
tant plasma effects are the thermal masses ~ gT acquired by the hard thermal particles [59l - [6l| . So a simpler (though 
less accurate) treatment for the regulator is to insert the Debye mass nip = (47ra ;s ) 1 / 2 T to the gluon propagator such 
that in the center-of-mass (CM) frame, 



\M 12 ^ 34 \ 2 CM a (127ra s ) 2 7V 2 
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where qr is the transverse component of q with respect to pi. It can be shown easily that Eqs. (|9]) and (fTTj) coincide 
in the CM frame in vacuum. This treatment was used in Refs. [3(3, [33|,[5(|. We will show both results using HTL and 
rriD for the mass of the internal gluon propagators for comparison. 

For the 23 process, because the matrix element is already quite complicated, we will just take as the regulator 
for internal gluons and estimate the errors. In the X^=iP« = convention, one can easily show that an internal 
gluon will have a momentum of ±(j»i +Pj) rather than ±(pi — pj). Therefore, the gluon propagator factors (ij) in the 
denominator of Eq. ([7]) , should be modified to 



1 



(ij) = oiiPi+PjY 



Pi ■ Pj 



2m 2 



(12) 



where we use m g to denote the external gluon mass. Then one applies Eq. ([8]) for the Boltzmann equation. In the 
numerator, the (ij) 4 combination is set by T and is 0(T S ). So we can still apply the substitution of Eq. (fT2"|) . even 
if the (ij) factors might not have the inverse propagator form. The error is ~ ^£)(*j) 3 = 0(a s T s ), which is higher 
order in a s . 

It is instructive to show that Eqs. (|7I8[) and (TT2|) give the correct soft bremsstrahlung limit. Using the light-cone 
variable 



P = (p + ,P ,PJ 



= (PO +P3,P0 -P3,Pl,P2) , 



(13) 



we can rewrite one momentum configuration in the CM frame in terms of p,p' , q and k: p\ — p, pi — p' , P3 = p + q — k, 
Pi = p' — q and p§ — k, with 

P = (Vi,m^,0,0), 
p' = (m 2 g /V^,V^,0,0) , 

q = (q + ,q^,q±) ■ 



(14) 



The on-shell condition p\ = p\ = m 2 „ yields 



-q±/vs, 

k]_ + yq\ - 2yk ± • q_L + (1 - y + y 2 )? 



Taking the large s limit, then the y — > 0, we obtain 
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where the permutation is over all final state gluon configurations. We see that Eq. (|16j)r educes to the Gunion-Bertsch 



formula [62j after taking mo,! 



0. A similar derivation of can be found in Ref. [63l. |64| . 



III. BEYOND VARIATION - SOLVING FOR r, SYSTEMATICALLY 

Following the standard procedure, the shear viscosity is related to B(p) by, 

N g f d 3 p 
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Bij (p) satisfies the constraint derived from the linearized Boltzmann equation, 
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x [By (p) + By(p 3 ) - By fa) - By(p 2 )] 

I dr 12;34p (i + ft q ){i + / 2 e9 )/ 3 e 7f (i + /p 9 )' 1 

4iv ff J1234 

X [By-(p) + By-(fc 4 ) + By(fc 3 ) - By' (fc 2 ) - Btffa)] 



E p 

® N 9 J 12:11 



rfr 123;4p (i + + / 2 e? )(i + / 3 e9 )/r(i + z; 9 )- 1 

x [By(p) + By(fc 4 ) - B«(fc 3 ) - Stf (*a) - Sii(fci)]. (18) 

However, solving By using this equation is technically challenging. It is easier to perform a projection (or convolution) 
to the above equation, then solve for the less restricted By. We will discuss the procedure below. 
By plugging Eq. (p~8|) into Eq. ([T7|) , we obtain 77 in a bilinear form of By , 



1 /" f/ 3 A- 
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x(i+/r 9 Ki+/2 9 )/3 e vr/r 

X [By(feg) + By(fc 4 ) + By(fc 3 ) - By (fc 2 ) - By (fc x )] 2 . (19) 

Then one can solve for By which equates the right hand sides of Eqs. (jTTJ) and (TIT)]) . This is nothing but a projection 
of Eq. (p~5|) . The resulting solution is not unique because only the projected equation but not the equation itself is 
satisfied. However, it is proven [28|, HH that the true solution of By, i.e. the solution satisfying Eq. (|T8|) . would give 
the maximum value of rj. Thus, solving for r\ becomes a variational problem. 

Recently, an algorithm is developed to find the true solution of By systematically Thus, this approach is 

no more variational but systematic. Here we outline the procedure. First, expanding B(p) using a specific set of 
orthogonal polynomials 0, EH : 

ax 

B{p) = (E p /TyJ^b r B^(E p /T), (20) 

r=0 

with y a constant chosen to be 1 in this case. The dimensionless polynomial B^ satisfying the orthonormal condition 

j0E~ f ^ q{1 + fp q )\P\ 2 ( E p/ T ) VB{r HE P /T)B^(E p /T) = T% s . (21) 
Then Eq. (|19[) can be written in a compact form, 

?"max 

n = (B \F\ B) = £ b r b s (b« |F| B^) , (22) 
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while Eq. (fT7|) gives 
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with 



L( o) = ^ B( o) I + f ;«)\ P f(E p /T)y. (24) 
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From Eqs. ([22]) and (l23l) . we can find bo by solving the equation 

r max 

L^S r0 = J2b s (B^\F\B^), (25) 



and then determine 77 from Eq. (f2"3")). [Note that there could be more than one solution satisfying Eqs. (|2"2"j) and 
but they all g ive the same 77.] 

In Ref. |66[, it is proven that this procedure gives a monotonically increasing value of 77 with increasing r max . Thus, 
one can systematically approach the true value of 77 by adding more terms in the expansion of Eq. (|2"U)) . We find good 
convergence in this algorithm. From r max = 1 to 2, 77 changes by less than 2% for a s < 0.3. Better convergence is 
found for smaller a,. 



IV. NUMERICAL RESULTS 
A. Leading-Log result 

The leading order [0(aj 2 )] contribution to 77 only involves the 22 and 23 processes [Hj]. The 22 collision rate is 
larger than 23 by a (lna s ) factor. In the leading-log (LL) approximation, one just needs to focus on the small qr 
contribution from the 22 process. Furthermore, it was shown in [671 . |68| that using the HTL regulator @ gives the 
same LL viscosity to that using the mjj regulator (JTTJ) . Thus, after performing the small qr expansion to Eq. (|19l) . 
we obtain 

77 LL ~ 27.1 T * (26) 

9 M\/g) 



which coincides with that of [27j to significant digits shown above. Using the entropy density for non-interacting 
gluons, s = Ng^-T 3 (for m g — 0), we obtain 

s - g'Hi/gy [ > 

This will be used to check our numerical result later. 



B. 7722 — Shear viscosity with the 22 process only 

To study the effect of the HTL regulator, 7722 (i.e. 77 with the 22 process only) with the HTL and tod for the 
internal gluon masses, respectively, are shown in Fig. [TJ The LL result till and AMY's 7722 [denoted as t]22(amy)] 

[27I [28( 1 are also shown. The external gluon mass m g , used in kinematics and in fp q such that E p = ^Jp 2 + m 2 , is a 

higher order effect in 7722- Changing m g from to mo = \f2rrioo yields an O^ri^/T 2 ) = 0(a s ) variation to 7722- This 
is confirmed numerically in the left panel of Fig. [TJ It is a good check to our numerical calculation that t\22(htl) an d 
V22(md) both converge to till in small a s , and r\22{HTL) agrees well with T)22(amy) when m g — is used to conform 
with the AMY result. 

About the HTL effect, r\22(HTL)l r \22(MD) is quite close to unity at a s = 10~ 6 , see the right panel of Fig. [JJ This 
ratio gets smaller at larger a s and reaches 0.65 at a s — 0.1 with little m g dependence (each T]22(HTL)/ r )22(MD) is 
evaluated with the same rn g ). This means the error is ~ 30% in the shear viscosity at a s = 0.1 if we use as the 
regulator for the gluon propagator instead of the HTL propagator. 



C. 7722+23 — Shear viscosity with the 22 and 23 processes 

In our full calculation, we use the HTL propagator for the 22 process. However, for technical reasons, we use the 
internal gluon mass nin for the 23 process. More specifically, we use matrix elements of Eqs. (|7 l9ll2p . E p — ^Jp 2 + m 2 
in kinematics and fp q for external gluons, and m g = rrioo. If the external gluons are massless but the internal gluons 
are massive, then the l/[(pi +P2) 2 ~ ntf)] factor could diverge. Using m g = m^, each term in Eq. ([7]) is non-negative. 
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FIG. 1: 7722 over the entropy density (left panel) and 7722 over rj 2 2(MD) (right panel) in various treatments. 'LL' is the leading 
log result of Eq. J26}. 'HTL' is the result using the full HTL matrix element of Eq. ([9]). 'MD' is the result using m_o as the 
regulator as in Eq. flT]). 'AMY' is AMY's result. 
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FIG. 2: Left panel: 1722/3 and 7722+23/5 for various cases. Right panel: the ratio of our result to AMY's. '22(HTL)+23' denotes 
7722+23 where m g — m^, the full HTL matrix element ([STJ is used for the 22 process, and the matrix elements of Eqs. (f7l l9H2[) 
are used for the 23 process. l 22(AMY)+23(AMY)' denotes AMY's result for 7722+23 (we have used the same s in 7722+23/s as 
ours with m g — m^). The range of the 'recommended value' of 7722+23 is bounded by 77- (lower bound) and 77+ (upper bound). 




In AMY and XG, external gluon masses were not included (m g = 0). This divergence was avoided by keeping only the 
most singular matrix elements in the small k±, q± limit, and taking the collinear approximation (AMY) or regulating 
the gluon bremsstrahlung infrared divergence by the Landau-Pomeranchuk-Migdal (LPM) effect (XG) which will be 
discussed in Section HV Dl 

We show 7722/s and 7722+23/s in the left panel of Fig. 2, where the HTL propagator is used for 22 and the external 
gluon mass m g = m^. We also show AMY's result for 7722+23/s for comparison. In the left panel, we see that our 
result for 7722+23/s deviates from XG's significantly, e.g. at a s — 0.1/0.01/0.001 our result gives 2.3/103/7020 but XG 
give about 0.45/20/1200 (read off from Ref. [3(|) respectively, i.e. our result is about 5 to 6 times as large as XG's. In 
Ref. [3l| they improve their calculation by using the Kubo relation and give larger values 0.795/60 at a s = 0.1/0.01, 
so our result is about 2 to 3 times theirs. In the right panel of Fig. 2, we see that the ratio of our result to AMY's 
approaches unity at a s < 10 -4 and ~ 0.8 at a s 0.004. The deviation in moderate a s is partly due to the finite 
angle, non-colinear 3-body configurations in the 23 process described by the full matrix element ([7]) and partly due 
to the gluon mass. We have also included a theoretical error band for 7722+23 which will be discussed in Section HVDI 

The effect of the 23 process can be seen more clearly in the ratio 7722+23/7722 shown in Fig. [3J We have plotted 
f]ii{HTL)+Til V22(htl) together with AMY's and XG's result for comparison. Our result shows that the 22 process 
dominates at small a s . When a s increases, 7722+23/^22 decreases and the central value reaches the minimum of ~ 0.6 
(which means the 23 collision rate is ~ 60% of the 22 one) at a s ~ 0.1 and then increases again for a s > 0.1. Thus, 
our result shows: (1) the 22 process is more important than 23 when a s < 0.1; (2) the finite angle contributions are 
important at intermediate a s (a s ~ 0.01-0.1). We see that AMY's result which employs the collinear approximation 
for the 1 «-> 2 process (corresponding to our 23 process), gives 7722+23/7722 close to unity. This implies their 23 collisions 
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FIG. 3: 7722+23 / '7722 in various cases. 'HTL' and uses HTL gluon propagators for the 22 process. The range of the 'recommended 
value' of 7722+23 is bounded by rj- (lower bound) and 77+ (upper bound). The external gluon mass is set to m^. AMY's and 
XG's 7)22+23 /f?22 are also shown. XG's result is taken from Ref. [3l| . 




10.6 



0.4 



0.2 



□ XG 

AMY 
— upper bound 
_22(HTL) + 23, m =rrv 

lower bound 



10" 



I 

10" : 



10' 



10" 




is just a small purturbation to the 22 collisions. XG's result in Ref. [30| which employs the soft gluon bremsstrahlung 
approximation, however, gives 7722+23 /V22 — [0-11, 0.161 around 1/8. This implies their 23 collision rate is about 7 
times the 22 one. Their improved result in Ref. [31| gives 7722+23/7722 — [0.1,0.3] corresponding to the 23 collision 
rate that is about 2 to 9 times (on average 5 times) the 22 one. Our result takes neither of the approximations, lies 
between AMY's and XG's result, but it agrees with AMY's result qualitatively but not with XG's. 



D. Error estimation 



Our T]22(ht x)+23 / ' s is tabulated in Table HI The error assignment is based on the the following error analyses for 
7722+23 : 

(a) HTL corrections for the 23 process. From our 7722 error analysis, we assign a ~ 30% error at a s =0.1 to the 23 
contribution for not implementing the HTL approach to the 23 collisions. The error will be smaller at smaller a s if 
the scaling for 7722 holds also for the error. Since the HTL effect tends to reduce the magnetic screening effect which 
lowers the IR cut-off and enhances the 23 collision rate, the HTL correction tends to reduce 7722+23- 

(b) LPM effect. We will try to estimate the error from neglecting the LPM effect. An intuitive explanation of this 
effect was given in Ref. [69| : for the soft bremsstrahlung gluon with transverse momentum fc-r, the mother gluon 
has a transverse momentum uncertainty ~ kx and a size uncertainty ~ 1/kr- It takes the bremsstrahlung gluon the 
formation time t ~ 1/ (fcyw^) ~ E^/kT 2 to fly far enough from the mother gluon to be resolved as a radiation. But if 
the formation time is longer than the mean free path l m f p « 0(a~ l ), then the radiation is incomplete and it would 
be resolved as gg — > gg instead of gg — > ggg. Thus, the resolution scale is set by t < l m f p . This yields an IR cut-off 
k\ > Ek/lmfp ~ 0(a s ) on the phase space [zjj. Thus, the LPM effect reduces the 23 collision rate and will increase 
7722+23- Our previous calculation using the Gunion-Bertsch formula shows that implementing the niD regulator gives 
a very close result to the LPM effect [33|]. Thus, we will estimate the size of the LPM effect by increasing the external 
gluon mass m g from 771,00 to mjj. 

(c) Higher order effect. The higher order effect is parametrically suppressed by 0(a s ), but the size is unknown. 
Computing this effect requires a treatment beyond the Boltzmann equation [38| and the inclusion of the 33 and 24 
processes. We just estimate the effect to be a s times the leading order which is ~ 10% at a s — 0.1. 

Combining the above analyses, we consider errors from (a) to (c). To compute a recommended range for 77/s, we 
will work with the -R22 and i?23 collision rates denned as 

= V22(HTL), 

(R22+R23) 1 = V22(HTL)+23- (28) 

Using HTL instead of mo in for the gluon propagator enhances the 22 rate by a factor of 

s _ R22(HTL) = V22(MD) ^ 
R22(MD) V22(HTL) 
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We will assume that the same enhancement factor appears in 23 rate as well, such that 

^ ^ R23(HTL) 
R23(MD) 

On the other hand, the LPM effect is estimated to suppress the 23 rate by a factor of 



(30) 



= R23 (Wg = m D ) 
R23 (m g = moo)' 

Combining the estimated HTL and LPM corrections to the 23 rate, the 22+23 rate is likely to lie in the range 
[R22 + R23, R22 + JSR23], while the higher order effect gives ±a s (R22 + R23) corrections to the rate. Without further 
information, the errors are assumed to be Gaussian and uncorrelated, the total rate is 

R22 + 2 ^-R 2 3) ± ( 2 ^ L R2 3 ) ± a s (R 22 + R 23 ) , (32) 



2 

and the recommended upper (77+) and lower (t]-) range for 7722+23 are 



V± = i 1 (33) 



R-22 + ^R23 ) T \l ( ^i?23 ) 2 + a 2 s (R22 + R23) 



2 



The rj± values are shown in the right panel of Fig. [2] and in Fig. [3l 

Our final 77/s result is as presented in Table I. At a s = 0.1, r\j s ~ [2.08, 2.58], which is closer to 2.85 of AMY than 
to XG's 0.5 in Ref. [30| and 0.795 ± 0.025 in Ref. [p. At a s = 0.3, we havej//s ~ [0.43, 0.759], which is compatible 
to 0.6 of AMY but not to XG's 0.13 in Ref. [H and 0.166 ± 0.025 in Ref. |31|. 



TABLE I: r\j s values for our ri22(HTL)+2s/ s result and the range of our 'recommended values' bounded by r)-/s and rj+/s. 





V22(HTL) + 23 
s 


V- 

« 


« 


a s 


V22(H TL) + 23 


1- 


v+ 

s 


1.00E-6 


3.75E+9 


3.75E+9 


3.75E+9 


0.100 


2.30 


2.08 


2.58 


1.58E-6 


1.54E+9 


1.54E+9 


1.54E+9 


0.125 


1.67 


1.48 


1.91 


2.51E-6 


6.33E+8 


6.32E+8 


6.34E+8 


0.150 


1.30 


1.13 


1.52 


3.98E-6 


2.60E+8 


2.60E+8 


2.61E+8 


0.175 


1.06 


0.909 


1.27 


6.31E-6 


1.07E+8 


1.07E+8 


1.08E+8 


0.200 


0.888 


0.752 


1.09 


1.00E-5 


4.42E+7 


4.40E+7 


4.44E+7 


0.225 


0.765 


0.638 


0.971 


1.58E-5 


1.83E+7 


1.82E+7 


1.84E+7 


0.250 


0.672 


0.551 


0.880 


2.51E-5 


7.57E+6 


7.52E+6 


7.61E+6 


0.275 


0.600 


0.484 


0.811 


3.98E-5 


3.14E+6 


3.12E+6 


3.16E+6 


0.300 


0.542 


0.430 


0.759 


6.31E-5 


1.30E+6 


1.29E+6 


1.31E+6 


0.325 


0.495 


0.386 


0.720 


1.00E-4 


5.43E+5 


5.38E+5 


5.47E+5 


0.350 


0.456 


0.349 


0.691 


1.58E-4 


2.26E+5 


2.24E+5 


2.28E+5 


0.375 


0.423 


0.318 


0.669 


2.51E-4 


9.47E+4 


9.38E+4 


9.55E+4 


0.400 


0.395 


0.291 


0.654 


3.98E-4 


3.97E+4 


3.93E+4 


4.01E+4 


0.425 


0.371 


0.268 


0.644 


6.31E-4 


1.67E+4 


1.65E+4 


1.68E+4 


0.450 


0.350 


0.248 


0.639 


1.00E-3 


7.02E+3 


6.94E+3 


7.11E+3 


0.475 


0.332 


0.230 


0.638 


1.58E-3 


2.97E+3 


2.93E+3 


3.01E+3 


0.500 


0.316 


0.214 


0.642 


2.51E-3 


1.27E+3 


1.25E+3 


1.29E+3 


0.525 


0.302 


0.200 


0.649 


3.98E-3 


542. 


532. 


553. 


0.550 


0.289 


0.187 


0.659 


6.31E-3 


234. 


229. 


240. 


0.575 


0.278 


0.175 


0.673 


1.00E-2 


103. 


99.9 


106. 


0.600 


0.267 


0.164 


0.691 


1.58E-2 


45.6 


44.1 


47.3 


0.625 


0.258 


0.153 


0.711 


2.51E-2 


20.7 


19.8 


21.7 


0.650 


0.250 


0.144 


0.736 


3.98E-2 


9.63 


9.08 


10.2 


0.675 


0.242 


0.135 


0.764 


6.31E-2 


4.62 


4.28 


5.01 


0.700 


0.235 


0.127 


0.795 
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E. When does the a s perturbation break down? 

We have only carried out the leading order [0(aj 2 )] rj/s in the a s expansion. Without computing the higher order 
contribution, it is hard to tell at what value of a s the perturbation starts to break down. In the above section, we 
have naively assumed the higher order contribution to be the leading order times a s , i.e. the expansion breaks down 
at a s ~ 1. However, explicit computations of thermal dynamical quantities and transport coefficients showed that 
the break down might happen at smaller a s (59l - [6l| (the screening mass computation breaks down at a s ~ 0.1 (oOj 
and the heavy quark diffusion constant computation breaks down at a s ~ 0.01 |6l|). 

Looking more closely to our leading order rj/s shown in the left panel of Fig. 2, the a s dependence of rj/s changes 
qualitatively at a s ~ 0.1. This could be a sign that higher order a s dependence has become as important as the 
leading order one. Thus, the higher order corrections could be bigger than our previous estimation, and our result 
might be only reliable when a s < 0.1. 

Having said that, it is interesting that our rj/s bends slightly upward at a s > 0.1 as if rj/s is trying to avoid 
going below the conjectured 1 /4tt bound. Several models have proposed to described the microscopic picture in the 
non-perturbative region |7lU74j[. And a similar result to ours is obtained in a recent calculation [75j based on one 
kind of simplification of the 23 matrix element [63[ . 

Our result implies that the proposed perfect fluid limit rj/s ~ l/(47r) cannot be achieved by perturbative QCD 
alone. 



V. CONCLUSIONS 



We have calculated the leading order [0(aj 2 )] rj/s of a gluon plasma in perturbation using the kinetic theory. The 
leading order contribution only involves the 22 and 23 processes. The leading order contribution only involves the 
22 and 23 processes. The Hard-Thermal-Loop (HTL) propagator has been used for the 22 matrix element, while the 
exact matrix element in vacuum is supplemented by the Debye mass mu for gluon propagators for the 23 process. 
Also, the asymptotic mass m.^ = mo/V^ is used for the external gluon mass in the kinetic theory as well. The errors 
from not implementing HTL and the Landau-Pomeranchuk-Migdal effect in the 23 process, and from the uncalculated 
higher order [(^(aj 1 )] corrections, have been estimated. 

Our result for rj/s lies between that of Arnold, Moore and Yaffe (AMY) and Xu and Greiner (XG). Our result shows 
that although the finite angle contributions are important at intermediate a s (a s ~ 0.01-0.1), the 22 process is still 
more important than 23 when a s < 0.1. This is in qualitative agreement with AMY's result. We find no indication 
that the proposed perfect fluid limit rj/s ~ l/(4ir) can be achieved by perturbative QCD alone. 

Acknowledgement: We thank G. Moore for providing us their tables of rj/s. JWC and QW thank KITPC, Beijing, 
for hospitality. JWC is supported by the NSC and NCTS of ROC. QW is supported in part by the '100 talents' project 
of Chinese Academy of Sciences and by the National Natural Science Foundation of China under grant 10735040. HD 
is supported in part by the Natural Science Foundation of the Shandong province under grant ZR2010AQ008. JD is 
supported in part by the Innovation Foundation of Shandong University under grant 2010GN031. 
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